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Fourier-transfrom Ghost Imaging based on Compressive Sampling algorithm 
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A special algorithm for the Fourier-transform Ghost Imaging (GI) scheme is discussed based 
on the Compressive Sampling (CS) theory. The CS algorithm could also be used for the Fourier 
spectrum reconstruction of pure phase object by setting a proper sensing matrix. This could find 
its application in diffraction imaging of X-ray, neutron and electron with higher efficiency and 
resolution. Experiment results are also presented to prove the feasibility. 

PACS numbers: 42.52.Ar, 42.50.Dv, 42.30.Wb,42.25.Kb 



The diffraction imaging of X-ray, neutron and electron 
provide significant methods to reveal microstructure p]- 
0] . Nevertheless, the traditional diffraction imaging by X- 
ray, neutron and electron in thermal state could only be 
used in periodic structure imaging. As for the aperiodic 
structure at the order of the source wavelength, accord- 
ing to the traditional wave theory, the coherent radiation 
source, such as Free Electron Laser (FEL), is essential for 
the imaging. However, since neutrons and electrons are 
fermions, their coherent sources with high brightness are 
in principle unavailable; although it's possible to obtain 
a coherent sources of photons with high brightness, the 
high requirement for brightness could only be met on the 
Synchrotron Radiation Facilities. 

It has been proved that two-photon correlation Ghost 
Imaging (GI) could realize the diffraction imaging by tak- 
ing advantages of thermal source |8[, however, it is also 
limited by the long acquisition time, as the correlation 
theory calls for mass samples to guarantee the ensemble 
average. Actually, lots of methods were introduced to 
improve its Convergence [9], however, as long as it still 
takes the correlation algorithm, the concept of ensemble 
makes it less effective. 

Recently the Compressive Sensing (CS) algorithm at- 
tracts more and more attention because of its extraor- 
dinary effect of reducing the samples [Iol-[l5|. GI and 
CS have similarities in extracting imaging: both are pre- 
pared with random sensing signal to "express" the imag- 
ing and a "point" detector to collect the result of expres- 
sion. But they also have intrinsic difference: GI is based 
on accurate "point by point" measure model [l6|, [I?}, 
while CS theory has proved that "global random" mea- 
sure model exhibits higher efficiency in imaging extrac- 
tion [10l-[l5j. Therefore, by combing GI and CS, it's pos- 
sible to develop a brand-new imaging model with higher 
imaging efficiency and resolution. 

At present, CS theory has already been introduced into 
GI scheme for real-space imaging [15j. Since CS algorithm 



could be generalized to complex field, it in principle could 
also be used to Fourier GI scheme, where most spectrum 
information relies on the phase part of the field. 

In view of these, we firstly propose a combination of the 
Fourier GI and CS in this letter. A recovery algorithm 
is investigated for the Fourier- Transform Ghost Imaging 
with Compressive Sampling (GICS) with experiment re- 
sults. This algorithm could not only improve the imaging 
efficiency to reduce the radiation damage of the sample, 
but may also improve the spectrum resolution compared 
to the Correlation Ghost Imaging (CGI). Therefore GICS 
may provide a brand-new diffraction microscopic imaging 
technic of X-ray, neutron and electron for the aperiodic 
structure. 

According to the CS theory, the algorithm is based on 
a corresponding relation between the imaging informa- 
tion X and the detect signal Y through a proper sensing 
matrix A EMI: 



AX = Y 



(1) 



If we sample K times, A should beaKxiV matrix, and 
Y should be a known K-element vector. Then by solving 
a convex optimization program, it gives an optimal result 
of X with N-element in the same expressing space as 

In our case, the lensless Fourier-transform GI scheme is 
shown in Fig 1 [8] . The thermal field from the source S is 
divided by the beam splitter (BS) into test and reference 
arms. To perform the Fourier-transform GI, there should 
be d\ = du2\ + <^22 for the scheme. Based on this scheme, 
we could express the intensity I w {r<i) on the test detector 
D2 as: 



Iw(r 2 ) oc J dxdx'E*(x)E(x')e(x)t(x') (2) 
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Where the integration of x and x r is over the object 
plane. By comparing with the CS theory, the object in- 
formation is located on the right side of (j2]), thus I w (r2) 
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could be considered as the known detect signal Y. Since 
our imaging goal is the Fourier-transform of t(x), to pro- 
cess the CS algorithm, it's necessary to establish a corre- 
sponding relation like (pQ) , and the sensing matrix A must 
be related to the signal I r (r\) from the reference detector 
Di, because the spectrum used to be expected to show 
up on the array detector D\(r\). 



FIG. 1: The lensless Fourier-transform GI scheme. Beam 
splitter (BS) is introduced to divide the field into two arms: 
in the test arm, field propagates freely cfei to an object and 
then c/22 to array detector D2; in the reference arm, field 
propagates dl to another array detector D±.di = d<z\ + d^i- 



If we consider the array of D\ to be large enough, then 
([2]) could be rewritten approximately as 



I w (r 2 ) oc J dndriJS7*(ri)JS7(ri) J dxdx' exp j-^p [{? ~ ri) 2 - (x f - r[) 2 ] J (3) 
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Compare (|3]) with ([T]) , the sensing matrix A could then 
be expressed as J E*(ri)^(r / 1 ) exp j— ^l\d 22 ( r i — r i 2 )}- 

Worth to note, unlike the sensing matrix A introduced 
in the bucket real-space GI scheme [15], the matrix A here 
is a K x TV 2 one rather than a K x N one. However since 
the reference detector Dl could only record N intensities, 
the N 2 — N interference terms of A when T\ ^ r[ are in 
general unobtainable . 

Generally, the only way available to get the exact field 
distribution on D\ in GI scheme is by calculating for 
the controllable source scheme or by introducing another 
known reference field E' for the homodyne detection. 
However, both are unavailable for the diffraction imaging 
of X-ray, neutron and electron. 

Notice if we consider each row of A to be a N x N 
matrix A sn 5 for each sample, the final imaging spectrum 
result obtained from CGI, |T(/)|, should lie just on the 
diagonal line of the CS solution X (also a TV x N ma- 
trix), corresponding to the known intensity distribution 



of Di on the diagonal line of A su b. According to the 
CS algorithm, the relative distribution of \T (f)\ could 
be expressed through the known I r {r\) of K sampling, 
while the existence of other non-diagonal elements just 
make the relation (|3]) hold. Therefore, if we only want 
to get |T(/)|, the accurate position of the non-diagonal 
elements, whose phases distribute randomly, become less 
important. 

Concretely speaking, if we properly conjecture the 
phase distribution on Di to be $(7*1) for each sample 
based on the concrete scheme, then each row of A could 
be written as 

A(r u r[) = ^// r (n)/ r (ri) exp {-i [$ (n) - $ (ri)]}4) 

XeXP {-A^^- ri2 )} 

Although the randomly conjectured phase dis- 
tribution [3>(ri) — $(7*1)] deviates from the true 



FIG. 2: Transmission function and spectrum intensity \T (f)\ 
of the pure phase five-slit object, (a) Transmission function 
of the object, the slit wide a is 150 /xm.(b) Spectrum intensity 
\T (/) | of the object observed from a lens of 5 cm focal length. 




value, the intensity distribution yjl r (ri)I r (r[) is ran- 
dom too, which makes the non-diagonal elements 
^/l r (ri)I r (r[) exp {— i [<£ (n) - $ (r[)]} in whole ran- 
domly distributed. When the number of conjectured val- 
ues is large enough, the values, though random, tend to 
cover all the true value of non-diagonal elements and es- 
tablish A su b. Since we don't care about T*(/i)T(/{), 
then follow the normal CS algorithm we could extract 
\T (f)\ along the diagonal line. 

Obviously, ((3]) is not a strict equation when the con- 
jectured values aren't enough, just as a non-ideal bucket 
detector for the real-space imaging [15j, whose influence 
to the sensing matrix is ignored. However, since we al- 
ways keep the true intensities I r (r\) on the diagonal line 
of A sn 5, this approximation has more influence on the 
efficiency to extract the information rather than destroy- 
ing the expressing basis for \T (f)\. Also, there's little 
relation in applications with strict equation for the CS 
algorithm to perform, thus it makes little sense to pursue 
the strict condition for CS algorithm. On the contrary, 
it's proper to develop the CS algorithm adaptive to rela- 
tions with different approximations. And still, GICS, as 
a brand new imaging algorithm instead of CGI, brings 
improvement of the low efficiency and limit resolution to 
the traditional CGI. 

Our experiment scheme is the same as shown in Fig.l, 
with source size a set to be 3 mm and wavelength is 
A = 0.532/im, c?2i = 20cm, ^22 = 5cm, d\ = ^21+^22, Di 
is an array detector and D2 is a point detector. The ob- 
ject is a one dimension pure phase five-slit, whose trans- 
mission function and spectrum intensity |T (/)| is shown 
in FIG. 2. The laser, CGI and GICS spectrum results 
of \T (f)\ are shown in FIG. 3 for comparison, where the 
GICS algorithm is based on the Spectral Projected Gra- 
dient for LI minimization (SPGL1) to solve the Basis 
pursuit (BP) problem in (pp) and we both take 100 off- 
diagonal elements of A sn 5 conjectured in FIG. 3(c) and 
(e). Obviously, from the same sample data GICS (shown 
in FIG. 3 (c) and (e)) obtains a much better result than 
traditional CGI result (shown in FIG. 3 (b) and (d)), 
which means improved extract efficiency. 

Besides the higher extract efficiency, GICS could also 
introduce super-resolution beyond CGI. This is of great 



FIG. 3: The experiment reconstruction results of spectrum 
intensity \T (f)\ comparison among laser, CGI and GICS. 
(a) Reconstruction results by laser illumination; (b) Recon- 
struction results by traditional CGI, the number of sample 
K = 100; (c) Reconstruction results by GICS (SPGL1), the 
number of sample K — 100 and number of iterations I = 30; 
(d) Reconstruction results by traditional CGI, K — 1000; (e) 
Reconstruction results by GICS (SPGL1), K = 1000 and 
7 = 30; 



importance to the reconstruction of object, especially for 
the Phase retrieval of nonperiodic objects [1[, where over- 
sampling method is available. It can be proved that the 
resolution of lensless Fourier-transform CGI is limited to 
the transversal coherent length on the reference detector 
D2. However, by taking advantage of GICS the limita- 
tion could be broken. 

The experiment for super-resolution is shown in FIG. 4. 
The scheme parameter is the same as FIG. 3, only with 
the object changed to a double-slit. The slit width is 
200/im, and the slit separation is 600/im. The transversal 
coherent length on D2 is calculated by l c = \di/ a ~ 
44/im, while the peak separation of the spectrum \T (f)\ 
is 32/im. Therefore the traditional CGI(shown in FIG. 4 
(b)) could not resolve the spectrum, and GICS (shown 
in FIG. 4 (c)) could realize super-resolution. The GICS 
result in FIG. 4 (c) is also based on BP problem in SPGL1 
but without any off-diagonal element of A sn 5 conjectured. 
This is also a proof for the role of the off-diagonal element 
discussed before, since it's still possible to extract \T (f)\ 
only based on the on-diagonal elements or the intensities 
on D2. In this experiment, both FIG. 4 (b) and (c) have 
been performed through spatial averaging technique. 

The improvement of extract efficiency and resolution 
introduced by GICS is also tenable for complicated spec- 
trum. FIG. 5 shows the GICS spectrum result of a com- 
plicated object "zhong" ring, compared with the results 
of laser illumination and CGI. The experiment scheme is 
also the same as in FIG. 3, and the diameter of the object 
ring is 0.8mm. Still the GICS result in FIG. 5 (d) is based 
on BP problem in SPGL1 without off-diagonal elements 
conjectured, and both CGI and GICS results (shown in 
FIG. 5 (b) and (c)) have been performed through spatial 
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FIG. 4: Spectrum super-resolution result of GICS. (a) Spec- 
trum result by laser illumination; (b) Spectrum result by 
traditional CGI with sample number to be K — 1000; (c) 
Spectrum result by GICS (BP problem in SPGL1), where 
K = 1000 and / = 500. 




(a) (b) (c) (d) 



FIG. 5: Spectrum result of a complicated object "zhong" ring, 
(a) Amplitude distribution of the complicated object "zhong" 
ring; (b) Spectrum result by laser illumination; (d) Spectrum 
result by CCS, where the sample number is K — 4000; (d) 
Spectrum result by GICS (BP problem in SPGL1), where 
sample number is K — 4000 and iteration number is / = 100. 



averaging technique. 

In conclusion, A new algorithm based on the lensless 
Fourier-transform GICS is presented in this letter. This 
scheme opens up a new way for the diffraction imaging 
of X-ray, neutron and electron with higher efficiency and 
super-resolution. The GICS, however, is much more than 
a algorithm revolution. From the point of microscopic 
essence, quantum fluctuations of field are generally de- 
scribed via high order correlation, which is based on ac- 
curate classical "point by point" measurement: the tra- 
ditional imaging technic is first order correlation (single- 
photon) system, while CGI is a second order correlation 
(two-photon) system. CS, on the contrary, stems from a 
"global random" sensing measurement. By introducing 
CS algorithm into GI, it also initiates a new discussion 
on its microscopic quantum mechanism. To find a rea- 
sonable explanation for the microscopic essence of GICS 
is a challenging problem, which will help to develop new 
imaging model with higher efficiency. 
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